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Abstract. The thermodynamics properties of globally neutral dense stellar 
matter are analyzed both in terms of mean field instabilities and structures be- 
yond the mean field. The mean field response to finite wavelenght fluctuations 
is calculated with the realistic Sly230a effective interaction. A Monte Carlo 
simulation of a schematic lattice Hamiltonian shows the importance of calcu- 
lations beyond the mean field to calculate the phase diagram of stellar matter. 
The analogies and differences respect to the thermodynamics of nuclear matter 
and finite nuclei are stressed. 
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1. Introduction 

Supernovae explosions, powered by the release of gravitational energy of a massive 
star which has exhausted its nuclear fuel, can lead to the formation of a most 
interesting dense stellar object: a neutron star[[rj. Due to the lack of observational 
data, the composition and structure of a neutron star is still highly hypothetical[|21. 
In the outer part of the star, the stellar crust extending over about one kilometer, the 
matter density is expected to be comparable to normal nuclear matter density, and 
the star can be modelized as essentially composed of neutrons, protons, electrons 
and neutrinos in thermal and chemical equilibrium. In a few minutes after its birth, 
the proto-neutron star formed at a temperature of « lO^^i^ becomes transparent 
to neutrinos and cools via neutrino emission to temperatures which are small on 
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the nuclear scale. The cooling process occurs via heat conduction and convection 
through the envelope to the surface on a time scale too short for the system to 
be in global thermal equilibrium, however local thermal equilibrium should be well 
verified during the whole evolution of the proto-neutron star, and beta-equilibrium 
is also often assumed[0. Crust matter is therefore very similar to nuclear matter, 
which is known to exhibit a complex phase diagram including first and second order 
phase transitions E] ■ 

The analogy between stellar matter and nuclear matter hides however an im- 
portant difference. If nuclear matter is by definition neutral, only global charge 
neutrality is required by thermodynamic stability for stellar matter. An important 
consequence of the charge neutrality constraint is that the canonical free energy 
density / is defined only for = pp - = 0. Hence f {T, pn, pp, pe) = f (T, Pn,p) 
with p {pp + Pe)/2 and the chemical potential pc associated to pc can not be 
defined since the free energy is not differentiable in the total-charge direction [ 
This constraint affects the thermodynamics directly since it changes the number 
of degrees of freedom of the thermodynamic potentials [ 6 . These considerations 
are expecially relevant when phase transitions are concerned. For exemple, since 
Pc is not a thermodynamic variable, the coexistence condition between two phases 
A and B p^ — p^ does not imply that each of the chemical potentials p^ and pp 
are identical in the two phases . The difference in chemical potentials of charged 
particles is counterbalanced by the Coulomb force: as some electrons move from 
one phase to the other driven by the chemical-potential difference, the Coulomb 
force reacts forbidding a macroscopic charge to appear. 

2. Mean Field approximation 

To illustrate the consequence of charge neutrality for the phase diagram of stellar 
matter, we present mean field calculations [C] using the Sly230a effective interaction, 
which has been optimized to describe exotic nuclei and pure neutron matter [1^. This 
force has been already applied to neutron stars crust in ref. [Ej. 

2.1. Homogeneous matter 

The matter instability to density fluctuations can be spotted looking at the curva- 
ture matrix [[Tn| 



where the free energy density is /(T, Pn,p) = fN{T, pn,Pp = p) + fe{T, Pe = p), and 
we have introduced the chemical potentials /i„ = dfN /dpn , Mp = (^In /dpp and pe = 
dfe/dpe- The additional term — dpe/dpe in the matrix modifies the stability 
conditions with respect to to the nuclear matter part, i.e. to the curvature of /at. 
Due to the small electron mass, i.e. high Fermi energy, the electron fluid is highly 
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Fig. 1. Minimum curvature of the mean field free energy matrix as a function of 
isospin asymmetry for normal nuclear matter(NM) and for stellar matter (NMe) 
at a temperature T = IQMeV calculated with the Sly230a effective interaction. 
Dashed line: stellar matter at zero temperature. 

incompressible leading to a quenching of the instability: the instability conditions 
DetC < or trC < are more difiicuh to fulfiU. The result for the Sly230a 
interaction is shown in Figure ^ the lowest eigenvalue of the curvature matrix is 
always positive independent of the temperature and proton fraction, meaning that 
the spinodal zone is suppressed in stellar matter [Q. 

2.2. Finite wavelenght fluctuations 

If we expect stellar matter to be stable at all temperatures respect to global density 
fluctuations, it is well known that the matter ground state in the inner crust should 
correspond to clusterized solid configurations composed of finite nuclei on a Wigner 
lattice[E]- To access these inhomogeneous configurations, stellar matter may be 
instable respect to finite wavelenght density fluctuations, 6pq = Aqt^^'"^ + A*e~'^^''^, 
with q = n,p,e. The mean field response to such a fluctuation induces two ex- 
tra terms in the curvature matrix eq.(|T)l[ IIT)]. The gradient term of the Skyrme 
fonctional[''5' produces a term oc fc^ coupling proton and neutron densities, which 
tends to suppress high k fluctuations; the direct Coulomb interaction produces a 
term oc fc~^ coupling proton and electron densities, which tends to suppress low k 
fluctuations. The resulting effect is the appearence of a flnite k interval where one 
eigenvalue turns negative, i.e. matter becomes unstable. An exemple is shown in 
Figure 12 12- We can see that an important instability region exists even at high 
temperatures and almost independent of the effective interaction employed, sug- 
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on the minimum eigenvaj^iaj^^s^ stellar matter curvature matrix at a fixedjV;^.jYp"^j 

of temperature, density ana proton fraction. Dashed line: value for homogeneous 
nuclear matter. Right side: instability region of homogeneous nuclear matter and 
inhomogeneous matter with electrons with two different effective interactions. 



gesting that clusterized configurations may be important also in the proto-neutron 
star evolution. 



2.3. The constraint of P -equilibrium 

If the cooling time scale is sufficiently slow, the proto-neutron star will also be 
subject to the constraint of /^-equilibrium, -I- /Lt*°* = /Lt^°* -I- /Lt^, with = 
IJ'q+rriql^^. This is certainly the case in the final step of the star evolution, when the 
temperature is low enough for the matter to be completely transparent to neutrino 
emission, /i^ = 0, which becomes then the most effective cooling process of the star[ 
Chemical equilibrium is generally assumed in the thermodynamic studies of 
stellar matter [EjEI and is often assimilated to charge neutrality. However we would 
like to stress that chemical equilibrium, as well as any other other transformation 
like constant proton fraction or constant temperature, is just a restriction of the 
accessible states and does not affect the thermodynamic properties, which are state 
functions [I13j. 

The chemical equilibrium constraint defines an accessible region in the phase 
diagram shown for a given temperature in Figure |3J 7 . The lower limit is given 
by the equation fi^ — while lepton number conservation defines an upper limit. 
We can see in Figure |31 that, depending on fii,, i.e. on the neutrino opacity, the 
instability region can be crossed under the constraint of chemical equilibrium. 
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potential (right side) pl^e ^^j^^^jixed temperature T — IQMeV. White rp^i(jj|y|j 

regions accessible to the constraint of /3-equilibrium. Some transformations corre- 
sponding to different neutrino chemical potentials are also given. 

3. Coulomb frustration beyond the mean field 

In the previous section we have seen that hot stellar matter can show an instability 
respect to finite wavelenght density fluctuations, eventually leading to the forma- 
tion of clusterized configurations. This instability is due to the interplay between 
the surface and coulomb interaction terms. This effect is a specific application of 
the generic physical concept of frustration. Frustration occurs in condensed matter 
physics whenever interactions with opposite signs act on a comparable lenght scales; 
applications range from magnetic systems to liquid crystals, from spin glasses to 
protein folding[ll4|. A well known application of frustration in nuclear physics is 
given by the multifragmentation observed in violent ion collisions. Concerning com- 
pact stars, it is recognized since a long time that in the clusterized state generated 
by frustration the structures may abandon spherical shapes and organize according 
to more complex topologies (pasta phases) [E]- The possible survival of these 
structures at high temperature is discussed in the recent literature[ ll6[ I17| . 

3.1. A statistical treatment of frustration 

A statistical tool to deal with frustrated systems is given by the multi-canonical 
ensemble[El- The two energy components En, Ec are treated as two independent 
observables associated to two Lagrange multipliers (3ntPc- A generalized grand 
potential is defined by 



(2) 
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Fig. 4. Event distribution dt '1^14,5 MM rfadd^iie'iife^;' Coulomb energy versus total 
energy plane for a finite nucleus of atomic number Z ^ 23 (left part) and Z = 82 
(right part) in the multi-canonical ensemble. Contour plots: uncharged system 
/3c = 0; contour lines: charged system Pc = Pn with different values of /Jat. Full 
line: spinodal curve. Insert: projection over the total energy axis for the uncharged 
(bimodal distribution) and the charged (monomodal distribution) case. 

where e"' is the fugacity of particle type q and the multi-canonical partition 
sum reads 



^;j„,/3c(^) = / dENdEcW(EN,Ec,N)exp {~(3nEn ~ (3cEc) ■ (3) 



If Ec represents the Coulomb energy and _Bjv the nuclear term, the choice [3c = 
(3n gives the usual (grand)canonical thermodynamics for charged systems, Pc — 
leads to the uncharged thermodynamics, while all intermediate values < Pc < Pn 
correspond to interpolating ensembles, or equivalently to physical systems with an 
effective charge (qeff/qo)'^ = Pc/Pn- The multi(grand)canonical ensemble allows 
to construct a single unified phase diagram for neutral and charged matter, and is 
therefore an ideal statistical tool to make some connections between nuclear and 
stellar matter. 

3.2. Application to multifragmentation 

As an example let us consider the phenomenon of multifragmentation as modelized 
in the Statistical Multifragmentation Model[ T8 . The event distribution in the 
multi-canonical ensemble eq.© is represented in Figure 01 for a value Pn corre- 
sponding to the transition temperature for the uncharged case. The first order 
phase transition gives rise to a two peaked distribution in the (multi)canonical en- 
semble, representing the two coexisting phases. The direction separating the two 
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peaks can be taken as a definition of the order parameter [(TQi . From the left part of 
Fig^lwe can see that the effect of the Coulomb interaction is a rotation of the order 
parameter which becomes almost orthogonal to the total energy direction, leading 
to a reduction of the latent heat and a shrinking of the coexistence zone. When the 
Coulomb effect becomes important (right part) , the spinodal region is not explored 
by the charged system and the phase transition becomes a cross-over. 

3.3. An Ising based model of stellar matter 




Fig. 5. Upper right: typical low temperature event distribution of the long ranged 
globally neutral Ising model. The proton-proton, proton-electron, and total dis- 
tributions are given. Left: coexistence region for the uncharged and charged case 
obtained from the multi-grand-canonical particle number distributions. Lower right: 
phase diagram. The constant proton fraction x = 1/3 line is also shown. 

To get a qualitative information of the effect of frustration on dense globally 
neutral stellar matter we have introduced a schematic but exactly solvable Ising 
model with long Coulomb-like and short nuclear- like range interactions [I2(J| 

<ij> i^j i^j 

where each site of a three dimensional lattice ol N = sites is characterized by 
an occupation number = 0, 1 and an effective charge qi = Ui — X^jLi "^j 1-^ ■ The 
effective charge represents the proton distribution screened by an uniform electron 
background [ 12] ■ To accelerate thermodynamic convergence, the finite lattice is 
repeated in all three directions of space a large number of times 
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Nn N 
n— 1 i,j 

with I = i + nL, J = j + nL. 




Fig. 6. Event distribution of the long ranged globally neutral Ising model in the 
Coulomb energy versus nuclear energy plane for different values of the Lagrange 
multipliers and different lattice sizes. 

The resulting phase diagram is shown in Figure [S] We can see that, at vari- 
ance with the case of multifragmcntation of finite charged nuclei, the effect of the 
Coulomb interaction is an expansion of the coexistence zone. This result can be 
easily understood inspecting the event distributions in the multicanonical ensemble 
in Fig|3[|2n|- We can see that the proton-proton and proton-electron energy dis- 
tributions are opposite in sign and similar in shape; in fact they would be exactly 
identical in the mean field approximation, meaning that the only effect of the charge 
at the mean field level is the addition of the electron free energy as we have discussed 
m section El The sum of the two contributions in this exactly solved model gives a 
two peaked distribution indicating coexistence between a high nuclear energy liquid 
and a low nuclear energy gas. The electron screening effect minimizes the Coulomb 
energy of the dense phase, making it accessible even at high temperature at variance 
with the finite nucleus case. This result implies that pasta phases may be relevant 
in a wide range of thermodynamic conditions . 

Another important effect of the Coulomb interaction in stellar matter is to 
suppress the critical point of neutral nuclear matter and the related phenomenon 
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of critical opalescence. The Coulomb energy density can be expressed as 



a f cTcjr, r') 
2V J |r-r' 



drdr' 27ra / ac{r)rdr 



(6) 



V 



where a{r,r') — {5pc{v)5pc{v')) is the charge density fluctuation. At the critical 
point o'c(r) is expected to scale as adr) oc r~^~^ , where is a critical exponent 
which turns out to be close to zero in most physical systems. A critical point 
would then correspond to a divergent Coulomb energy. As a consequence, the 
phase transition in stellar matter ends at a first-order point. This can be clearly 
seen in our numerical simulation as shown in Figure 1^ I2(J|. The critical point of 
the uncharged system is characterized by a diverging coulomb energy (upper right), 
while for the charged system at the end point of the coexistence zone (lower right) 
the fluctuation does not increase with the lattice size. 

Such an effect has been already observed in Ising models with long-range frus- 
trating interactions, where the coexistence region was seen to end at a flrst-order 
point [ I14| . Since the Fourier transform of the correlation function gives the en- 
hancement of the in medium scattering cross section respect to its free value, the 
important physical implication of this result is that we expect hot stellar matter 
to show a small opacity to neutrino scattering [ EI j which may have important 
consequences on the cooling dynamics [I22|. 

4. Conclusions 

In this contribution we have estabilished some connections between exotic phases 
in stellar matter and nuclear multifragmcntation. Both phenomena are triggered 
by an instability to finite wavclcnght density fluctuations and can be understood in 
terms of frustration between the nuclear and the Coulomb interaction. 

If the first order phase transition of normal nuclear matter is suppressed in stel- 
lar matter due to the high electron incompressibility, a (second order [ 6 ) transition 
is however expected between single-phase and coexistence configurations. Large 
correlated structures are therefore predicted, even if at the ending point of the 
coexistence region the correlation lenght is not expected to diverge. 

Since the Coulomb energy is minimized in homogeneous partitions correspond- 
ing to pure phases, and it is maximal in clusterizcd partitions, a higher tempera- 
ture is needed to reach the mixed partitions with respect to the uncharged system. 
Therefore, the mixed-phase phenomenology may be relevant for the proto-neutron- 
star structure in a wider temperature range than usually expected [151. Such an 
expansion of the coexistence region is at variance with the mean-field expectations . 
This contrast stresses the importance of realistic calculations beyond the mean-field 
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